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A framework which combines Green's function (GF) methods and techniques from the theory 
of stochastic processes is proposed for tackhng nonhnear evolution problems. The framework, 
established by a series of easy-to-derive equivalences between Green's function and stochastic 
representative solutions of linear drift-diffusion problems, provides a flexible structure within 
which nonlinear evolution problems can be analyzed and physically probed. As a preliminary 
test bed, two canonical, nonlinear evolution problems - Burgers' equation and the nonlinear 
Schrodinger's equation - are first treated. In the first case, the framework provides a rigorous, 
probabilistic derivation of the well known Cole-Hopf ansatz. Likewise, in the second, the ma- 
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chinery allows systematic recovery of a known soliton solution. The framework is then applied 
to a fairly extensive exploration of physical features underlying evolution of randomly stretched 
and advected Burger's vortex sheets. Here, the governing vorticity equation corresponds to the 
Fokker-Planck equation of an Ornstein-Uhlenbeck process, a correspondence that motivates an 
investigation of sub-sheet vorticity evolution and organization. Under the assumption that weak 
hydrodynamic fluctuations organize disordered, near-molecular-scale, sub-sheet vorticity, it is 
shown that these modes consist of two weakly damped counter-propagating cross-sheet acoustic 
modes, a diffusive cross-sheet shear mode, and a diffusive cross-sheet entropy mode. Once a 
consistent picture of in-sheet vorticity evolution is established, a number of analytical results, 
describing the motion and spread of single, multiple, and continuous sets of Burger's vortex 
sheets, evolving within deterministic and random strain rate fields, under both viscous and in- 
viscid conditions, are obtained. In order to promote application to other nonlinear problems, 
a tutorial development of the framework is presented. Likewise, time-incremental solution ap- 
proaches and construction of approximate, though otherwise difficult-to-obtain backward-time 
GF's (useful in solution of forward-time evolution problems) are discussed. 

Keywords: Green's function methods, stochastic methods, advection-diffusion problems. Burger's 
equation, nonlinear Schrodinger equation, Burger's vortex, random strain field, Feynman-Kac solu- 
tion of vorticity transport, Cole-Hopf derivation, Ornstein-Uhlenbeck model of vorticity transport, 
hydrodynamic organization of vorticity 



1. Introduction 

Parabolic evolution equations of the form 

A{x, T, 'n)'n^ + B{x, T, 'n)-Vr] + C{x, r, r])V^ri = F(x, r, r?) (1) 

lie at the heart of a far-reaching set of physical theories and models. Here, rj, the variable of 
interest, and A, B, C, and F can represent scalar, vector, or higher order tensor quantities, and 
X and T correspond respectively to space- and time-like variables. A limited list of examples 
include the Navier-Stokes equations (see, e.g., [Ij), Schrodinger's equation [2], advection-diffusion 
equations [3], reaction-diffusion equations [3], equations governing nonlinear pattern formation 
[5] , a diverse set of wave equations [HIT] , Chapman-Kolmogorov equations governing evolution of 
Markov processes [8], and particle and continuum versions of mass, charge, and linear and angular 
momentum conservation PlIlUj. 

Finding general approaches for treating nonlinear versions of ([T]) drives a large field of re- 
search |1H - I17J . Numerical methods are typically favored in applied problems, while intense interest 
attaches to development of analytical techniques since these can deepen physical and mathematical 
insight [TTl[lll[16l[l7], and somewhat secondarily, allow, e.g., code benchmarking and algorithm 
development [7l[T8]. 

Green's function (GF) methods, of course, find wide application to linear ([1]), where the method's 
versatility derives from several features: i) GF solution structure exemplifies simplicity, providing 
a transparent description of a linear system's response to often complicated space- and history- 



dependent boundary and initial conditions, and internal forcing; ii) GF's employ, via delta func- 
tions, space- and time- localized, but physically-lumped i.e., black-box, descriptions of typically 
ill-understood system interactions with these forcing agents; and iii) GF's can often be ascribed 
intuitive interpretations [51 ll9l - [25] . promoting physical understanding. 

A number of papers have used Green's function approaches to tackle nonlinear evolution prob- 
lems |1 H I18 1 [2^. However, with the exception of a limited collection of systematic techniques, e.g.. 
Green's element methods [18] and the inverse scattering transform pLlj, work in this area remains 
in an early state of development. 

A well-known, though little-used stochastic representative solution for the non-homogeneous 
Cauchy problem, a linear drift-diffusion embodiment of ([T]), was presented in 1931 by Kolmogorov 
|25p26|. In analogy with the Green function's circumvention of unresolved small-scale forcing, 
the stochastic solution rests on Wiener's [8] physically lumped description of unresolved random 
forcing. The latter feature, paralleling the role played by delta functions in GF models, promotes 
tremendous versatility in both the reach and interpretation of stochastic process models. 

This paper pursues four objectives: 
A. Given the linear nonhomogeneous form of equation dTJ, subject to nonhomogeneous initial and 
boundary conditions, Kolmogorov's stochastic representative solution can be expressed [251126]. 
Our first objective centers on highlighting a series of equalities that exist between expectations 
appearing in Kolmogorov's solution and corresponding terms in the Green's function solution. This 
demonstration is significant since it establishes a simple bridge between two broad fields. Green's 
function methods and the theory of stochastic processes. 



Importantly, recognition of this connection establishes a structure in which stochastic and 
Green's function methods can be applied in concert to a range of deterministic and random evo- 
lution problems of the form in ([T]). We refer to this structure as the Green's function- stochastic 
methods framework (GFSM). 

B. Given this framework, the remainder of the paper focuses on two principal questions: 

a) Can the framework be systematically applied to solve nonlinear versions of equation ([1])? 

b) Can the framework facilitate physical and mathematical exploration of problems characterized 
by some element of randomness? 

With regard to the first question, since most nonlinear problems require some form of numerical 
attack, our second objective focuses on set-up of time-incremental solutions. The main ideas include: 
i) derivation of exact, space-dependent Green's functions, valid over arbitrary, though small time 
increments, ii) identification of forward and backward time evolution problems and associated 
adjoint problems, and iii) as illustrated in Test Case 1, use of approximate forward-time Green's 
functions as surrogates for hard-to-compute backward time GF's. 

C. The third objective, also addressing question a), centers on testing the viability of time- 
incremental attacks on nonlinear versions of ([T]). Two nonlinear evolution equations having ex- 
act solutions. Burger's equation and the nonlinear Schrodinger equation, provide test beds. In 
both Test Cases, we start from a general, time-incremental Green's function solution, and identify 
strategies for obtaining known non-incremental solutions. 

D. The last objective, addressing question b), aims at illustrating how GF and stochastic process 
ideas can provide essential physical guidance in the analysis of nonlinear evolution problems. For 



this demonstration, and as presented in section 6, we study a canonical fluid mechanics problem 
which captures the ubiquitous, combined effects of vortex stretching, advection, and diffusion: 
evolution of Burger's vortex sheets. 

As described in section 6, much of the analysis pivots on gaining an understanding of near- 
molecular scale, sub-sheet vorticity transport and organization. Crucially, this question emerges 
when we observe that evolution of sheet-scale vorticity, initiated by a delta function initial condition, 
corresponds to evolution of the transition density for an Ornstein-Uhlenbeck process. 

In other words, in making this connection, we are presented with a radically alternative, prob- 
abilistic picture: development of any given Burger's vortex sheet, evolving within either a deter- 
ministic or random strain rate field, can be interpreted as the collective, stochastic evolution of a 
swarm of elemental vortex sheets (EVS). 

Thus, in order to identify a reasonable physical embodiment of EVS's, we are lead to investigate: 
i) the highly disorganized structure of sub-sheet, short-time-scale vorticity, and ii) the long BVS- 
time-scale organization of this vorticity. 

Mathematically, couching analysis of BVS evolution in terms of the stochastic evolution of el- 
emental vortex sheets likewise proves advantageous since it allows straightforward, probabilistic 
determination of time-dependent BVS mean position and spread. Similarly, in the case of con- 
tinuous initial vorticity distributions, a stochastic vantage point leads naturally to Feynman-Kac 
solutions for the resulting vorticity evolution. 

In closing the Introduction, and as an aid to navigating the paper, we highlight the paper's 
essential five-part structure: 



I) Mathematical details needed for application of the framework to linear and nonlinear evolution 
problems are given in sections 2 and 3. 

II) Testing use of time-incremental GF's for solving nonlinear problems is described in sections 
4 and 5. 

III) Application of Green's function and stochastic process ideas as physical and mathematical 
probes is, as mentioned, illustrated in section 6. The paper's final three parts correspond 
to three distinct elements comprising the illustration. Thus, set-up and calculation of single 
sheet GF's is carried out in sections 6.1 and 6.2. 

IV) Sections 6.3 shows that evolution of individual Burger's sheets can be interpreted as the 
stochastic evolution of a swarm of sub-sheet elemental vortex sheets, the latter representing a 
quasi-physical embodiment of an underlying Ornstein-Uhlenbeck process. Section 6.4 pursues 
a physical interpretation of this observation, and in the process, addresses the fundamental 
question of how highly disorganized, near-molecular-scale vorticity becomes organized on long 
BVS length and time scales. 

V) Finally, sections 6.5 through 6.7 investigate the evolution of single, multiple, and continuous 
collections of Burger's vortex sheets advected and stretched by random strain rate fields. 

2. Green's function-stochastic methods framework 

As mentioned, the GFSM framework can be viewed as a union of Green's function methods and 
the theory of stochastic processes, bridged by the equalities ([HI) through (fT6]) derived below. This 



section focuses on derivation of equations (|14p through ()16p . Section [3] and the Appendix describe 
time-incremental apphcation of the framework to nonhnear problems. 

Although it is likely that the approach used here - derivation of the Green's function solution 
followed by comparison with the backward-time stochastic solution - has been described elsewhere, 
we have not located such descriptions. We note that Friedman [2B], following Kolmogorov |25] . 
derived an essential result: the representative stochastic solution of a backward time, linear, non- 
homogeneous advection-diffusion problem. However, [25li26] did not connect the stochastic solution 
to the equivalent Green's function solution. 
Here, we use the following non-rigorous recipe: 

i) Define an appropriate forward time linear evolution problem. 

ii) Define the associated adjoint problem governing the Green's function, G. 

iii) Use the adjoint problem to derive the Green's function solution. 

iv) Write the linear evolution problem in backward time form and recognize that the result is a 
backward Kolmogorov (Fokker-Planck) equation; thus, express the solution as a representative 
stochastic solution [26]. 



v) Finally, equate corresponding terms in both solutions, i.e., equate terms involving boundary 
conditions, the initial condition, and the nonhomogeneous forcing term. 

Remark 1: It is important to note that for divergence- free drift fields, the validity of this ap- 
proach can be easily proven: under these circumstances, the adjoint problem governing G, and 
the backward time Fokker-Planck problem governing the transition density, p, are identical. Thus, 



the Green's function can be interpreted as a transition density, and hence the Green's function- 
stochastic solution equahties in (jl4|) - (J16p become identities. 

Remark 2: For apphcation of the GFSM to nonlinear evolution problems, the appropriate linear 
evolution equation simply corresponds to a linearized version of the original. As illustrated in 
Test Cases 1 and 2, linearization takes advantage of the fact that the solution is constructed time- 
incrementally. 

2.1. Evolution problem template 

The examples treated in this paper can be mapped in some fashion, e.g., time-incrementally for 
nonlinear problems, to the following linear drift-diffusion problem: 

Mr? = -/(x',i') onQ = Dx{0,t] (2) 

7/(x', 0) = (/)(x') on 6Qo = Dx {t' = 0} (3) 

77(x', t') = g{ji.', t') on 6Q = 6D x (0, t] (4) 

where the operator, M, is given by 
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., dr] 



Mr] = i^V^r] + h- V'r] - —^ (5) 

and where, for notational convenience, we use the backward drift b in place of the forward (actual) 
drift V (with b = — v). In addition, 0(x') is the initial condition, g{x',t') is a time-varying Dirichlet 



condition, the differential operators V ^ and V' denote derivatives taken witii respect to x', and D 
is the spatial solution domain. See the Appendix for a description of forward and backward time 
coordinates and forward and backward evolution equations. 

2.2. Representative stochastic solution 

The backward time stochastic solution of the evolution problem in ([2])- ([5]) within a space-time 
domain, Q = D x [s,T), subject to Dirichlet conditions, rj = g{x,s'), on the boundary 5Q of Q, 
and a final condition rj = (f>{x,s = T) on the final time-slice, D x {s = T}, can be expressed in 
representative form as 



r^{x,s) = E^^,[g{x{T),T)]+E^^,[(l){x{T))]+E^,s[f f{x{s),s)ds'] (6) 

where s = T — t, is backward time, t forward time, and T is the backward time instant at which 
the field r/(x, s) is known. We outline a rough proof of this below. Here, -Ex,s is the expectation 
associated with the stochastic process, x(s), sampling respectively, Dirichlet conditions on 5Q, 
the final time condition, (f){x,T) on D x {s = T}, and the forcing function / within Q. Thus, r 
represents the random time at which the process meets (and is absorbed on) the Dirichlet boundary. 

2.3. Bridge relations between Green's function and stochastic solutions 

In order to clearly highlight the origin of the bridge relations, we first outline well-known steps 
for generating the GF solution to ([2]) - ([S]), and then present a non-rigorous, though easy-to-grasp 
derivation of Kolmogorov's solution, equation ^. 
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Green's function solution 

The Green's function solution slightly extends the well-known derivation outlined, e.g., by Morse 
and Feshbach [28j and Barton [29j; the extension consists of incorporating non-zero drift, v = — b 
into the governing equation for r]. 
The derivation, stated in tutorial fashion, is as follows: 

i) Given the forward time operator M defined in ([5]), find the corresponding adjoint operator 
M* (see, e.g., [30]), here given by |26J: 



M* = i^V'^ - b • V' - V' • b + — (7) 

ov 



ii) Define a Green's function, G(x,t|y, t'), satisfying 

M*G = -6{t - t')6ix - y) (8) 

for all X, y G D, and 0<t' <t. 
iii) Multiply (^ by G(x,t|x',i') and dSD by ri{x',t'), and form: 

/ / [G{Mr]) - r]{M*G)]dD'dt' = - I [ [Gf - i]5{x - x')6{t - t')]dD'dt' (9) 
Jo Jd Jo Jd 

which allows application of Green's theorem. Here, e/t << 1, and G = for t' > t. 
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iv) Carry out the integrations in ([9]) to obtain the solution for r/(x, t) : 

r?(x,i)= / I GfdD'dt'-u [ i rjV'G-n'dS'dt'+ f (l){x')G{x,t\x',0)dD'- [ i GJ-n'dS'dt' 
Jo Jd Jo Jsd Jd Jo Jsd 

(10) 
where J = vrj — u'V'rj is the total flux of rj (and where again b = — v). 

Stochastic solution 

This derivation, again presented in recipe fashion, is a wholly heuristic, non-rigorous version of that 
given by, e.g., Friedman [26]. A non-technical approach is preferred since it provides a simple path 
to the desired result, equation ([6]). 

i) Consider the differential, random, multidimensional displacement of the stochastic process, 
X, over the backward time interval ds' : 

dx{s') = h{x{s'), s)ds' + \/2^(iw(s') (11) 

where w is a multi-dimensional Wiener process. For any given realization of the random 
displacement, dx(s'), use a Taylor expansion about the solution point, (x, s), to compute the 
associated change in rj : 

dT]= ris'ds' + h-V'ri + uV^r] ds' (12) 

where, in anticipation of taking the expectation over the Wiener process, terms in dwi = 
■Wi{s + ds') — Wi{s), are set to and terms in dwidwj are expressed as ds'6ij. In addition, we 
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keep ds' small enough that quadratic and higher order terms in ds' can be neglected. Finally, 
since dwi is a zero- mean, gaussian random variable, expectations of terms involving dw^ 
are zero for n = 1, 2, 3, ... See, e.g., Gardiner [S| for further details. 

ii) Next, integrate p^ along the random path traced by the process x(s') as it progresses from 
s toward the backward time, T : 

ri{x,s) ='I]{x{t),t) - rjs'ds' + h-'V'ri + i/V'^'i] ds' (13) 

J s 

where T corresponds to the forward-time initial instant (say, t = 0), and where the desired 
term, r/(x, s), has been isolated. Again, in anticipation of taking expectations, terms involving 
dwi, dwidwj with i ^ j, and dw-^"' have been excluded. Here, the time r depends on where 
x(s') ends up: for those realizations that impact the hyper-surface, 5Q = 6Dx [0, T), enclosing 
the space-time solution domain, Q = D x [0,T), prior to reaching the final space-time slice, 
D X {s' = T}, T is the (random) time of impact. For those realizations that survive without 
impacting 6Q, t = T. 

iii) Next, replace the argument in the last integral with the right side of ([2|), take expecta- 
tions with respect to the process %, and use the fact that ii^x,s^(x (''"),''") = -E'x,s5 {x{t)jt) + 
E:x,,s4> ix{T)) , to obtain the final result, equation ([6]). 

In order to obtain the bridge relations, we simply recognize that in linear problems, t] at any 
point (x, s), as determined by the representative stochastic solution in ([6]) must correspond to 
r]{x,t), as computed via the Green's function solution in p^ . Thus, comparing similar terms in 
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([6]) and ([TO]) , noting the correspondence between the backward solution time, s, and the forward 
solution time, t, we obtain the following: 



E^,s g{x{r)) =-^ f i 5(x',t')V'G(x,t|x',t') • ndS'dt' 

Jo JSD 



(14) 



^x,. <t^{x{T)) = I (/>(x')G(x,t|x',0)(iZ)' (15) 

Jd 

E^^s r f{x{s'),s')ds'= [ /" /(x',t')G'(x,t|x',t')dx'dt' (16) 

Js Jo Jd 

where dH) has been used in the second term on the right in (jlOp . 

Note that an equivalent form of the the equality in (jlSp . valid for Cauchy problems in unbounded 
domains, is rigorously proven in Friedman [26]. Likewise, a relation equivalent to (J14p . and appro- 
priate to the Dirichlet problem Mrj = 0, is given by Schuss [27J. In addition, in the case where drift 
b = —V is everywhere zero, the solution in pup is identical to that given by Barton [2U]. Finally, 
stochastic representations are available for Neumann and mixed initial boundary value problems; 
see, e.g., [31]. 



3. Incremental solutions for nonlinear problems 

When confronting nonlinear and/or nonhomogeneous problems, use of a time-incremental attack 
is immediately suggested. The idea is simple - shrink the forward solution time interval. At' (or in 
backward evolution problems, the backward interval As') in order to: i) allow linearization of non- 
linear evolution problems, and ii) when necessary, allow nonhomogeneous sources, F(uj_|_i,x, Tj+i), 
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to be expressed as F(uj,x, Tj), where in either forward or backward time, Tj+i = Tj + At. 

In this section, and in anticipation of Test Cases 1 and 2, we write down incremental versions 
of the stochastic and Green's function solutions in (^ and pUj) . The Appendix describes essential 
requirements that incremental GF's must meet. 

3.1. Incremental Green's function and stochastic solutions (1-D case) 



Test Cases 1 and 2 are initiated from an incremental version of (jlOp . In both, we assume that 
the solution point (x, t) € Q is sufficiently removed from boundaries to allow neglect of boundary 
conditions. Thus, referring to (|1U|) . setting D = (— oo, oo), and dropping boundary terms, we arrive 
at 



ri{x,tj+i)=l G{x,tj+i\x',tj)rj{x',tj)dx' + / G{x,t\x',t')f{x',t')dx'dt' (17) 



where, to allow incorporation of linearized nonlinear terms as well as nonhomogeneous terms, the 
nonhomogeneous source, /, is included. Physically, (fT7|) gives the response, r]{x,tj+i), at tj+i due 
to both the 'initial condition', r]{x',tj) at tj, and the time-dependent forcing, f{x',t'), that takes 
place over Atj = tj^i — tj. 

For completeness, we also write the associated time-incremental version of the representative 
stochastic solution in ([6]): 

Vix,tj+i) = E,,t.^,rj{x{tj),tj) - E^^t,+, I ' f{xis'),s')ds' (18) 



15 



where, for notational consistency, we have stated the backward time coordinate, say Sk, in terms 
of its corresponding forward time coordinate, tj+i (thus, for example, Sk -^ ij+i and Sk+i -^ tj). 

Finally, note that incremental versions of the bridge relations in (|14| ) - (|16| ) can be easily written 
down by comparing like terms in, e.g., p7|) and (|18|) . 

4. Test Case 1: Solution of Burger's equation 

In brief overview, this example tests application of the GFSM framework to nonlinear drift-diffusion 
problems. As an aid to understanding the development, we note the following essential points: 

i) A minimal requirement for application of time-incremental GF's to nonlinear problems rests 
on derivation of an incremental GF. By properly limiting the time step size, the difficult-to- 
solve backward-form adjoint equation can be recast in approximate, soluble, forward form. 
Section 4.2 assumes a simple stochastic process viewpoint in order to derive generic time step 
constraints, equations (p7|) and ([28j) . that must be met when using this procedure. 



ii) Nonlinear problems having a space or time-dependent drift, u, lead, not unexpectedly, to 
incremental GF's containing the same. See ()3ip below. As will be shown, it is sometimes 
possible to transform the nonlinear GF to linear form by eliminating the a priori unknown 
u. Here, two steps are used: 

(a) Assume a transform, 0(n), having a specific parametric form in u, and require that (p be 
governed by a simple (linear, non-advective) diffusion equation. The form of cj) used is 
given by equation ([3T 
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(b) Insert the assumed (j){u) into the diffusion equation and compare the result (equation 
()38p ) with the original nonlinear governing equation ( ()19|) below); attempt to determine 
a detailed form for </>(u) by forcing the former equation to match the latter. 

Here, and in response to |38j who noted the arbitrary origin of the Cole-Hopf solution to Burger's 
equation, we find a suitable (j){u) and thus, a fundamental derivation of the Cole-Hopf solution. 
Burger's equation, 

Ut + UUx — I'Uxx = (19) 



which gained prominence as a simplified model of compressible irrotational flow [32], has since found 
application in a wide range of other problems, including shock dynamics [6j, nonlinear acoustics [50], 
magnetohydro dynamics [H], turbulence [12], traffic flow [l3], dynamics of dislocations, polymer 
chains, and vortex lines |44J , and formation of large scale cosmic structure [45j . Variants of Burger's 
equation also appear in models of flame front propagation and forced diffusion [46j. 
In the following, the source-free Kardar-Parisi-Zhang (KPZ) equation 



ht + -hl- vhxx = (20) 



appears as a useful analytical intermediary. In brief, this equation serves as a well-studied model 
of interfacial growth ^47j, where the local rate of change in interface height, /it(x, t), reflects surface 
smoothing, vhxx, due to, e.g., condensation and/or evaporation, combined with growth, h^., in the 
local surface-normal direction. Here, we only require the well-known transformation between the 
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KPZ and Burger's equations: 

u = h^ (21) 

(where in multiple dimensions u = V/i); see, e.g., [35]. In this case. Burger's equation describes 
the nonlinear evolution of the interfacial slope. 

4.1. Time step conditions allowing construction of approximate incremental G 
for backward-form adjoint problems 

The adjoint equation in G is, as shown below, of backward form, and thus proves difficult to solve. 
However, by restating the adjoint equation in forward time form and by properly choosing the 
time step size. As, we can obtain an approximate analytical solution. This section presents a 
straightforward argument for determining appropriate conditions on As. The conditions obtained 
are general, applicable to any evolution problem in which advection and diffusion are extant. 

To begin the incremental solution over t' G (tj,tj+i], replace both r]{x',t') and —b{x',t') = 
— b-Cx = u{x',t') in ©-(Is]) with u{x',t'), and set f = g = 0. The associated adjoint equation, from 
([7]) and dH), then assumes the form: 

Gf + u{x', t')G^' + Gu^'ix', t') + z^G^/^/ = (22) 

with 

G{x,tj+i\x',t') ^ 5{x' - x) as t' ^ tj+i (23) 

where again tj < t' < tj+i. 



18 



Restate (|22p (which is in backward form in the forward time t') in forward form, in the backward 
time, s' : 

Gs' + b{x', t')G^, + Gb^, {x', t') - uG,.^, = (24) 

with 

G(x, Sfc|x', s') — )■ 5(x' — x) as s' ^ s^ (25) 

where s^ < s' < Sfc+i) and where s^ and s^+i correspond respectively to tj+i and tj. Note, 
As = Sfc+i — Sfc and At = tj+i — tj. 

Next, determine conditions on As that, first, allow b{x',s') in (j24p to be approximated by 
b{x,Sk), and second, allow neglect of the term Gb^f, given these conditions, (|24p can be linearized 
and an analytical solution obtained. 

We proceed heuristically. First note that since G satisfies the same (Fokker-Planck) equation 
and initial condition as an incremental transition density function, p{x, Sk\x', s'), associated with a 
stochastic process 

dx{s') = bixis'), s'){s' - Sfc) + V2^[wis') - wisk)] (26) 

then for small As, the random walk swarms governed by (|26p and launched from the x — axis at 
s' = Sfc+i, toward the backward time slice s' = Sk (i.e., the current forward time, t) : i) will at 
Sfc_|_i, have an approximate mean launch position x{sk+i) ~ x + b{x,Sk)As, and ii) will have non- 
negligible probability of reaching the solution point (x, s^) only if they lie within an approximate 
distance vi^As of x(sfc+i). Hence, the size of the region over which both G = p and gradients in G 
are non-negligible, is on the order of \^uAs. 
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Denoting the length and time scales associated with the drift field b as Xs and tg, respectively, 
and noting that the respective scales of the four terms in ([M|) are 0(Gs/As), 0{bsGs/\^iy^s), 
0{bsGs/xs), and 0(Gs/As), where Gg and bg are the local Green's function and drift scales, then 
we find that the first condition on As, allowing neglect of the Gb^' in (|24|) . is 



ViyAs , , 
« 1 (27) 



A second condition. 



As , , 

— « 1 (28) 



allowing replacement of b{x'^s^) with b{x^Sk)^ follows by expanding b{x\s') about (x^Sk) • 



bix',s')=bix,s,) + 0[^^^,^] 

Xs i^s 

where x' — x = 0{y/iyAs) and s' — s = 0(As). 

Assuming that As is chosen so that (j27p and (I28p hold, then (j24l) can be solved via Fourier 
transform 

/•CO 

G{x,Sk\k,s')= e-2^^^^'G(3;,Sfc|x',s')dx' (29) 



oo 



with the result 



r.( I / /^ 1 i -[x'-{x + b{x,sk){s'-sk))]' 

G{x,Sk\x,s)= ^ ^ ^=exp ^ -— '— (30) 

y/ATru{s' - Sk) L 4z/(s' - Sk) 



This can be restated in terms of physical (forward) time variables using s' — s^ = ij+i — t' and 
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b{x,Sk) = —u{x,tj+i) = —u{x,tj) + 0{At), to arrive at 

nr , I ' ,^ 1 i -[x'-{x-u{x,t,){tj+i-t'))]^ ^ 
G[x,tj+i\x ,t ) = , exp — -— (31) 

Thus, using (j3ip in (J17p (with / again set equal to 0), we obtain the incremental Green's function 
solution to Burger's equation: 

u{x',tj) J 4i'7r{t j ^i —tj) exp — [{x' — Xj+i)'^/{4u{tj+i — tj))]dx' (32) 

-oo 

where Xj+1 = x{tj^i) = x—u{x,tj){tj+i—tj). Since the incremental Green's function, G{x,tj+i\x',t') 
corresponds to an incremental transition density, p{x,tj+i\x',t'), then this solution can also be in- 
terpreted probabilistically as an explicit version of the representative solution in (llSp . 

4.2. Derivation of the Cole-Hopf solution 

Although the presence of the term u{x,tj) in the exponential in (|3'2p does not prevent, e.g., 
numerically-based solutions, in order to both validate the time-incremental approach as well as 
explore potential approaches for analytically tackling other nonlinear evolution problems, we now 
focus on using ()32p to obtain the non-incremental Cole-Hopf solution. 

The path connecting (|32p to the Cole-Hopf solution is indicated by the form of the incremental 
solution. In particular, we recognize that finding a non-incremental solution requires, at minimum, 
elimination of the a priori unknown drift, u{x',t'), from the incremental Green's function in (|3ip ; 
elimination of the drift is required in order to stretch the time increment. At', to arbitrary lengths. 
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Thus, we introduce a surrogate, (p = (p{u), for u upon which we impose two key requirements. 
First, in order to take advantage of the machinery developed above, we require that (f>{u) also 
satisfies an incremental evolution equation of the same form as 



oo 



x,tj^i) = / (f){x' ,tj)G'''{x,tj-^-i\x',tj)dx' (33) 



where G^ is the incremental Green's function associated with the evolution of (p. Second, in order 
to eliminate u from the incremental Green's function, C^, we require that the evolution of be 
purely diffusive: 

(pf - i^<Px'x' = (34) 

so that the associated adjoint equation is 

Gf, + i/Gf,^, = -6{t' - t)5{x' - x) (35) 

The Green's function, C^, is easily obtained by setting u = in (|3ip : the detailed incremental 
solution in (j33p then follows: 

(j){x,tj+i) = cl){x',tj) JATTu{tj+i-tj) ex.p-[{x' -x)'^/{Av{tj+i-tj))]dx' (36) 

where now x = x. 

In order to proceed, we next recognize that by assuming a specific parametric form for 0(n), 
followed by introduction of this guessed form into (I34p . an evolution equation in u emerges. If the 
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latter can be transformed or forced to match the actual equation for u, (J19p . then we will have 
constructed a transform from the incremental to the non-incremental solution. Thus, after a few 
trials, we choose 

^{u) = exp(/(u)) (37) 



Inserting (|37j) into (|3l|) then yields 



f" 

f/ , J 1„ 2 



Uf - fu^'x' -v[f' + —]u^ = (38) 



where /' = df/du. 

No choice of f{u) allows (j38p to be placed in the form of Burger's equation (J19p : however, by 

redefining u as h, and choosing 

f" 1 

^' + 7 = -^ (39) 

we observe that (|38p transforms to the KPZ equation (j20p . Replacing u with /i in (|38p and differ- 
entiating the result with respect to x' then yields 



Vf - i^vx'x' - Ivvvx' \r + 77] - ^"^^ if + 7]' = ° (^°) 



where, for clarity, we express h^' as v, and where the last bracketed term is differentiated with 
respect to u. Thus, by (|39p . ()40p transforms to Burger's equation (J19p . and as noted, (j38p transforms 
to the KPZ equation i^O\i . 

Although a number of solutions for f{u) via (|39p are available, we choose /' = — l/(2z^) [where 



23 



it is understood that u in ([37|l now corresponds to h in (pOj) and w in (jl0|) corresponds to u in p^ ]. 
Thus, from dST]) 

</)(/!) = 0(x',O=exp(^) (41) 



or inverting, 



/i(x',t') = -2uln{<i){x' ,t')) 



(42) 



which corresponds to the Cole-Hopf transformation for the KPZ equation j48) . 

The corresponding Cole-Hopf transform for Burger's equation, which we denote as (piv), then 
follows from (|4ip via integration oi v = hx '■ 



{v) = 4>{x',t') =exp(- 



/o^ v{x",t')dx\ 
2^ ' 



(43) 



Returning to the incremental solution (|33p for (p (with (p is restated as (p), we insert (|43|) and 
follow Whitham [S] by isolating v{x,tj^i) : 



v{x,tj+i) = -2u 



x,t 



•j+i 



(44) 



or 



v{x,tj+i) = 2v- 



V47ri/(t,-+i - i,)] no</'(a;',ii)exp-[(x'-x)V(4i.(t,-+i -t,-))]2];g^^rf^' 



V47ri/(tj+i - tj) /!^(/.(x',t,)exp-[(a;'-x)V(4z.(tj+i -t,))]dx' 



(45) 



Finally, set tj+i = t, tj = 0, and after some manipulation, obtain the classic Cole-Hopf solution of 
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Burger's equation (see, e.g., [6] ): 






where R = H{x'; ,x,t) = J^ v{y, t = 0)dy + ^^ ^ ^ 



2t 



5. Test Case 2: Soliton solution of the nonlinear Schrodinger equa- 
tion 

As a second test, we consider solution of the one-dimensional nonlinear Schrodinger equation 

iVt' + r]x'x' + K|r/pr/ = (47) 

Like Burger's equation, (j47p represents a canonical nonlinear evolution equation [6j, capturing 
in this instance nonlinearly dispersive, weakly dissipative wave propagation. As with Burger's 
equation, ()47p and its variants appear in a wide range of contexts, including nonlinear optics [6], 
hydrodynamics [49j, and plasma physics |50j . 

Here, we wish to show how a well-known soliton solution to (j47p can be obtained using the 
machinery developed above. In contrast to the function transform approach used in Test Case 1, 
however, we make the jump from incremental to non-incremental solution via asymptotics; such 
approaches can be considered, e.g., when nonlinear terms in the governing evolution equation are, 
on the scales of interest, small. 
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To begin, we seek a traveling wave solution [6j of the form: 

r/(x', t') = expiigx' - ift')h{X') (48) 

where X' = x' — Ut' is a coordinate attached to the moving wave, U is the wave speed, and g and 
/ are constants. 

We focus on determining the wave envelope shape, h{X'), and simply note that g and / can be 
determined in terms of U using, e.g., Whitham's approach [6j. Thus, (j^Tj) is re-expressed in the 
wave-fixed coordinate system as 

ir]t' - Ur]x' + Vx'X' + K\r]\'^r] = (49) 

The corresponding Green's function, again valid over short time intervals, At', is governed by 

- iGt> - UGx' + Gx'X' + K\r]\^r] = -6{X - X')6{t - t') (50) 



In order to determine G, we take the spatial Fourier transform of (|50p . 



/oo 
exp{-2TTikX')G{x, t\X', t')dX' 
-oo 



to obtain: 



1 r°° 

G{X,t\X',t') = — exp[27rik{X' -X)]ex]i[[-2TTUk + iATr^k\t' -t)]dk (51) 

^'^ J -oo 
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The incremental solution over At' then follows from ([T 



+ IZo I-oo 2^e2-^fc(^^') exp [[-27TUk + iA7r^e]At'] h{X')e^^'-f^'+^''UkdX' (52) 

where t' = t' - t and AX' = X' - X. 

Focusing first on the first term on the right of (I52p and carrying out the time integral, we note that 
since {2tt)-^ j°<^^ ^2nikAX' ^ s{AX') = 5[X' - (x-Uot)], then this term becomes e^(»^-^*)/i3(X)AT, 
where, for clarity, we write Uq = —U, and where At = At'. 

Turning to the second term in (I52p and using similar steps, we arrive at 



/oo poo -1 p -. 

/ _e2^^^^^'/i(X') exp i [g{X' + [/„(* + At)) - f{t + At)] dX'dk = 
-oo J —oo ^'^ '- 

/oo P 

5{AX')h{X') exp i[g{X' + [/o(t + Ar))(iX'/i(j; - Uo{t + At)) exp [i[gx - f{t + At))] (53) 
-oo 

where AX' = X' — {x — Uo{t + At)), and where the term exjp[27rUok + i47r^A;^] At has been replaced 
by 1 + 0(At). 

The incremental solution (1521) thus assumes the form 



h{x - C/ot)e^(s^--^*) = KATh^ix - C/ot)e^(9^--^*) + h{x - Uo{t + AT))e*(»^--^*)e^^^^ (54) 



Canceling the common exponential term and expanding the last term above about X = x — Uot, 
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then yields 



h{X) = KArh^iX) + e'f^^ [h{X) - hxAX + ^-hxxAX^] + 0{5X^) (55) 



or equivalently, 



^^—h{X) = KArh^iX) + -hxxAX^ - hxAX (56) 



where AX = UoAt. 

In the last expression, we observe a separation of scales, i.e., the term involving hx is of 0(At), 
while all remaining terms are either of order Ar^ or O(kAt). This becomes clear by expressing 
(j56p in dimensionless form 

hj^j^AX^ + fih^ -pJi- h<^AX (57) 

where li = h/hs, X = X/xs, AX = AX/Axs, h = 2kAt/i^x^/Ax^, J3 = fdr'^xl/Axl, Axs = 
UqAt, and where hg and Xg denote the amplitude and axial length scale of the wave envelope, 
h{X). [Note, AX is 0(1).] 

Thus, expressing h for example, as 

h = ho + ehi + 0{e^) (58) 

where e = Axg/xg, and focusing on the case where k and /3 are both positive and 0(1), we obtain 
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a physically and mathematically consistent solution for h{X) : 



ho = (59) 

and 

h = ejy^ecH^^rx (60) 

which is a well-known soliton solution |6j of the cubic Schrodinger equation (|47p . 

As a closing summary to this and Test Case 1, we have shown that a systematic attack on 
nonlinear evolution problems can be initiated from an incremental Green's function solution. In 
most problems, one would typically proceed to a numerical time integration, using the incremental 
GF to construct the kernel. Here, for purposes of validation, analytical integration has been 
pursued. 

6. The GFSM as a physical probe: organization of near- molecular- 
scale vorticity in Burger's vortex sheets 

The last example studies physical features underlying evolution of single, multiple, and continuous 
sets of Burger's vortex sheets evolving within deterministic and random strain rate fields. The 
example is designed to illustrate application of Green's function and stochastic process ideas as 
probes for exploring linear and nonlinear evolution problems. Since the example is long, encom- 
passing three distinct elements, we expand upon the introductory overview. 

The generic problem of evolution of N^ Burger's vortex sheets (BVS), generated by the ap- 
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pearance of an arbitrary, spatially discrete, one-dimensional velocity field is first presented (section 
I6.ip . A discrete initial velocity condition is chosen since it leads to a physically and mathematically 
crucial delta function initial condition on the vorticity evolution problem. 

Mathematically, as detailed in section 6.2, this condition allows calculation of an essential single 
sheet Green's function; the importance of this GF emerges when it is reinterpreted as a transition 
density and applied to compute statistical behavior of single, multiple, and continuous sets of BVS's 
evolving within random strain rate fields (sections 6.5-6.7). 

Physically, and as mentioned, the delta function IC allows us to interpret the vorticity transport 
equation as a Fokker-Planck equation for an underlying OU process (section 6.3). Presuming that 
the OU process describes the stochastic dynamics of a sub-BVS physical entity, viz, elemental vortex 
sheets, we are led to investigate vorticity dynamics within individual Burger's sheets (section l6.4p . 
This examination leads to the following observations and results: 

1) On short acoustic time and near-molecular-length-scales, scaling shows that in-sheet vorticity 
is disordered, three-dimensional, and diffusive. 

2) Since vorticity, on the long BVS time-scale, becomes one-dimensional and highly organized, 
some organizing mechanism clearly operates over the longer time scale. 

3) Presuming that organization is effected by weak, in-sheet hydrodynamic modes, we investigate 
these modes using a simple analog: sub-sheet vorticity organization in unstrained planar 
vortex sheets. 

4) Finally, the modal analysis suggests that three hydrodynamic mechanisms underlie in-sheet 
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organization: 

(a) weakly-damped, cross-sheet acoustic modes, 

(b) a diffusive cross-sheet shear mode, and 

(c) a diffusive cross-sheet entropy mode. 

Burger's vortices were first studied by Burgers in 1948 [39j who showed that such structures are 
capable of dissipating turbulent kinetic energy through the combined action of viscous dissipation 
and vortex stretching. Soon after, Townsend [51] proposed that the fine structure of high Reynolds 
number turbulence, i.e., the structure extant on scales smaller than the viscous dissipation scale, is 
characterized by random distributions of Burger's line and sheet vortices. Since these early works. 
Burger's vortex lines and sheets have been studied both as a fundamental, analytically tractable 
model of the combined action of viscosity, advection, and stretching on vorticity, see, e.g., 
and as a putative fundamental structure in various turbulent flows 



6.1. Initial conditions and governing equation 

Attention is limited to the case where a series of parallel vortex sheets are formed at some instant, 
t = to, within a two dimensional potential flow, i.e., a strain rate fleld, given by u = [—kx,0, kz\. 
See figure 1. The formation of each sheet occurs due to the appearance of a spatially varying flow 
component in the y-direction, given by 

v{x, t = to) = 2, ^ViU{x — Xi) (61) 
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Figure 1: Continuous collection of Burger's vortex sheets evolving under influence of an inviscid 
strain rate field, Vx = —kx, Vz = kz. Vorticity is one-dimensional in the z— direction and is produced 
by the y — velocity component, Vy{x,t). 

where v = Vy, Avi is the incremental velocity at Xi, U{x — Xi) is the unit step function, and Nx is 
the number of increments. There is no limit on the size of individual velocity increments; as shown 
immediately below, these determine the strength of each associated vortex sheet. 
The vorticity equation assumes the linear 1-D form 



i^t — kxil,x = k^ + l^Oa 



(62) 



where Q is directed along the z — axis. The initial condition associated with (j6ip is 



N^ 



0(3;, t = to) = y^ Avi6{x — Xi) 



(63) 



i=l 



Due to the linearity of the governing equation ()62p . the last condition can be given a clear physical 
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N:, 



meaning: it corresponds to the superposition of N^, infinitesimally thin vortex sheets, {rijjj^ 



n{x,t = to) = ^^i{xi,to) (64) 



where 

^i{xi,to) = /^Vi5{x - Xi) (65) 

(and where units on 5{x — Xi) = -^ are length"^). 

The generic solution, describing the vorticity field produced by the A^^^ sheets, is of the form: 

n{x,t) = ^^i{x,t) (66) 

where the detailed form of Vti is given below. 

6.2. Derivation of the single sheet Green's function 

We derive the single sheet Green's function, G^*'(x, i|xj, to), by first normalizing ()62p with the initial 
sheet strength, /\vi , and then by noting that the resulting problem takes the form of that governing 
a forward time GF; see, e.g., |29j. Thus, we solve for r2j(a;, t)/Afj, and use 



G«(x,t|xo,to) = ^M^ (67) 



where Xq is the initial location of the infinitely thin vortex sheet. We take to = 0, temporarily 
drop superscripts and subscripts referring to sheet i, and refer to the vorticity field produced as an 
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individual Burger's vortex sheet. 

While Townsend's 1951 paper [51] is cited by Saffman [52j as the source for a solution to a scalar 
transport problem occurring in a time-varying strain field of the type considered here, Townsend's 
derivation, based on the method of characteristics (MOC), is apparently given elsewhere. Here, 
using a Fourier transform-MOC approach, representing a slightly generalized version of Gardiner's 
solution for a one-dimensional OU process evolving in a steady drift field [_8j, we first derive a 
solution appropriate to time-varying strain fields of the form: 

V = —k{t)xex + k{t)zez (68) 

Fourier transforming (|62p yields 

Clt + k{t)KVL^ + vK^Vi = (69) 



Placing this in characteristic form then yields: 



-— = -z/K^O on — = k{t)K (70) 

iJLL (JjL 



Integrating along characteristics from Kq = n{to) to k = K,{t) yields 



«;(i) = Koexp[/i(t)] (71) 
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where 

ft 



h{t) = / k{t')dt' (72) 

'to 



Next, substituting ([71]) in ([TUD gives 



il{K, t) = Cl{ko, to) exp [-i^kl I e^^^*''^dt'] (73) 

Jto 

while setting the initial time to = and using the initial condition 

n{x,Q) = /:^v5{x-Xo) (74) 

leads to 

Vl{ko, 0) = /^ve'""^'" = Av exp[me-''Wxo] (75) 

Finally, letting 

p(t) = e-2Mt) TeSMOdi' (76) 

Jo 

and taking the inverse transform of ()73p leads to the final single-sheet Green's function: 



Gfx, t\Xo, to = 0) = — = — , exp — (77) 

^ ' ^ Av ^J^^up{t) ^ ^Mt) ^ 



The single sheet GF shows that peak response to the initial delta function travels to x{t) = XqC"'*'*' 
and spreads as p{t). As a quick check, we note that in the case where k{t) is constant, (j77|) assumes 
the appropriate form [52','53]. 
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In order to circumvent destabilizing vortex compression, it appears that for most times, t, 
k{t) must remain positive. Although short periods of vortex compression, sandwiched between 
long periods of stretching, can likely be stably sustained, we make no attempt to address this 
question. Rather, we limit attention, in the deterministic case, to k{t) > 0, and in the random case, 
k{t) = ko + k'(t), to positive ko, with ko > \k'{t)\max, where ko and k'{t) are the non-random and 
random parts of k(t). 

6.3. Vortex sheet evolution as an Ornstein-Uhlenbeck process 

Focusing on the evolution of an individual BVS, say the i*^ sheet (having vorticity Qi{t)) in a 
collection of Nx sheets, again defining a normalized vorticity, fij = ^i/Avi, and again noting the 
initial condition (j65p . we observe that equation ()62p can also be interpreted as a Fokker-Planck 
equation governing the transition density, 

p»(x,t|x„t,) = ^M^ (78) 

Avi 

of an OU stochastic process X (O? where x evolves as: 

dx^^^ (t) = -k{t)x'^^ {t)dt + V2^dw^^ (t) (79) 

Importantly, this connection allows us to introduce a correspondence between the evolution 
of individual realizations of the stochastic process, X (t), and evolution of individual elemental 
vortex sheets. In this picture, the instantaneous BVS, here the i sheet, corresponds to a cloud 
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Figure 2: Decomposition of individual Burger's vortex sheets into elemental vortex sheets. 

of say N elemental vortex sheets; see figure 2. Although we suppose for the moment that EVS's 
are quasi-physical entities, analogous to say idealized, infinitely thin vortex sheets, we argue below 
that a physically reasonable embodiment of these can be defined; see section I6.4.2I 

Given a stochastic picture of BVS evolution, we can now easily write down physically meaningful 
expressions for the mean position of the vortex cloud (= i Burger's vortex sheet) < x (^) >x, 
and it's diffusive spread, i.e., its thickness, < {x (t)— < X (t) >x)'^ >x'- 



<X^\t)>, 



xp^'"' {x,t\xi,to)dx 



(80) 



< ix^'\t)- < x^'^it) >x? >x= / (x- < x^^^t) >.fv^'^{xAxi,to)dx (81) 

J —oo 

where expectations over the stochastic process, X (0; ^^^ written in the notationally convenient 
form, < • >a; . 
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When the strain rate k{t) is random (see sections 6.6 and 6.7), equations ([62]) . (f65|) . ([79 |) -([8T ]) 
apply to any reahzation of k{t). Thus, fohowing, e.g., [55H57] . the average sheet position and spread 
over random k will be computed as: 



( < x^'^(t) >x) = { I xp^'\x,t\xi,to)dx) (82) 

and 

< < ix^'Ht)- < X^'\t) >^f >.) = { I {x- < X^'Ht) >,fp^'\x,t\x,,to)dx) (83) 



where expectations taken with respect to k will be denoted as < • > . 

6.4. Elemental vortex sheets 

In order to develop a physically and mathematically consistent picture of EVS's, we limit attention 
to incompressible flows and consider vorticity transport on length and time scales that are small 
relative to those associated with BVS dynamics. 

In order to identify an appropriate time scale, we first appeal to the Biot-Savart law: 



1 f rxcj(y,t) , , ,, 

vs(x,t) = -- /^^dy (84) 



where r = x — y. This kinematic result, which gives the instantaneous solenoidal (incompressible) 
velocity field, Vg, in terms of an integral over the vorticity field, u, indicates that the local magnitude 
and direction of the former instantaneously senses and responds to the integrated effects of the 
latter. In reality, sensing and response are acoustically mediated. 
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Thus, we surmise that in order to expose sub-sheet dynamics, we should focus on processes tak- 
ing place on acoustic time scales. This choice, in fact, proves advantageous since the acoustic scale 
lies well separated from both the much longer BVS time scale, k~^, and the much shorter molec- 
ular collision scale, Tc, given below; as will be shown, substantial insight into sub-sheet rotational 
dynamics emerges on this intermediate scale. 

For the EVS length scale, we focus on lengths that are large relative to molecular scales, do, and 
small relative to the characteristic BVS thickness, di{t), where, depending on whether the strain 
rate field is deterministic or not, 6i{t) is given, respectively, by ([HT]) or ([83]) . 

Thus, as a means of gaining a conceptual foothold, it proves useful to focus on the acoustic 
time scale, near-molecular length scale rotational dynamics and evolution of in-sheet clumps, fiuid 
particles comprised of a fixed number of say Nd molecules. Limiting attention to liquids, and given 

1/3 

a characteristic (effective) molecular diameter do, the characteristic clump size is 6ci ~ Nj do- 

For a short period, Td, any given clump remains nominally intact and experiences an incremental 
rotation, the magnitude of which corresponds to the average rotation of all constituent molecules. 
Since clump dispersion occurs due to thermal motion of constituent molecules, Td ~ Sd/vT, where 
vt = \/3kT/m, the thermal speed, is on the order of the sound speed, Oq- 

1/3 

Since the ratio of Td to the collision time scale, Tc = 0{do/vT), is 0{Nj ), choosing Nd to 
be on the order of say, 10^, gives Td/Tc = 0(10). Hence, as depicted in figure 3, we envision that 
an initially smooth collection of clumps becomes 'bumpy', i.e., thermally roughened, over a time 
interval of O^Td); likewise, as noted in the caption, individual clumps remain largely intact. 

Indeed, the Td — scale inter-clump transfer of rotational momentum can be viewed as a conse- 
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Figure 3: Over a time interval on the order of Td, an initially smooth collection of neighboring 
clumps begins to thermally disperse. Since the clump is large relative to individual constituent 
molecules, dispersion on this time scale takes place predominantly via thermal roughening and 
inter-penetration of surface molecules (where the latter are depicted as surface bumps). Holes 
within, and penetrating jets into, clumps form as well, but the fraction of molecules having sufficient 
energy to form these over Td is small; thus, these are not depicted. 
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Figure 4: On acoustic time scales and on the elemental vortex length scale, vorticity corresponds to 
the instantaneous average rotation of N^i molecules comprising a fluid clump. On these short time 
and length scales, and under subsonic conditions, the vorticity field is three-dimensional, disordered, 
and strongly diffusive. 

quence of the cog-like action of clump-surface molecules interacting with one another. Figure 4 
provides a schematic representation of vortex sheet vorticity on the large BVS scale, the smaller 
elemental vortex scale, and on the clump-scale. 

Qualitative insight into clump-scale dynamics can be gained by assuming that a continuum 
description applies on scales of order 5ci and Td- Recognizing that clump-scale gradients in BVS-scale 
velocities and vorticities are small and nondimensionalizing the full vorticity transport equation, 
we obtain: 

Cbi - Re-^V^cD = -Ma-it • V • cD (85) 

where Re~ = y/{ao5ci) is a clump-scale Reynolds number. Ma = Ug/ao is a BVS-scale mach 
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number, Us = 0{kxs) is the BVS velocity scale, Xg is the BVS-scale sheet position, and where 
we have dropped a 0{5d/xs) term describing stretching of BVS-scale vorticity by the clump-scale 
velocity field. 

Note first that for liquids like water, and for clumps having 0(10'^) molecules. Re = 0(1). Under 
these circumstances, and under conditions where the BVS-scale mach number is small, (j85|) shows 
that, relative to diffusion, clump-scale stretching (again, dropped from (j85p ) and advection of vor- 
ticity are likewise small. This is consistent with a detailed analysis of organization of clump-scale 
vorticity by in-sheet hydrodynamic modes below (section 16. 4. ip : there, it is shown that diffusional 
smoothing via shear and entropy modes constitute two of three hydrodynamic organizing mecha- 
nisms. 

Thus, from (|85p . the following picture of clump-scale vorticity transport emerges: 

a) in low speed incompressible flows (Ma << 1), transport is three dimensional and diffusive, 

b) advection, important on BVS length and time scales, only emerges on the clump-scale when 
liquid bulk velocities are high, on the order of the sound speed, and 

c) stretching remains weak to the point of nonexistence. 

6.4.1. Organization of clump-scale vorticity: hydrodynamic modes 

Given that vorticity transport on BVS scales is highly organized and one dimensional, evolving 
on these scales under the combined action of advection, stretching, and diffusion, while sub-sheet 
(clump-scale) transport is highly disorganized, three dimensional, and diffusive, it becomes apparent 
that some mechanism, acting over time and length scales long relative to 5ci and t^i, organizes 
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disordered clump-scale vorticity. 

We assume that the organizing mechanism is associated with weak hydrodynamic modes su- 
perposed on the bulk, BVS-scale velocity field. In order to expose the role of these modes on 
sub-sheet organization, we focus on a simpler analog problem, organization in unstrained vortex 
sheets {k = 0). 

The analysis adapts Mountain's |58i well known approach for analyzing the small wave number, 
low frequency response of simple liquids to excitation via scattered inelastic light and neutron 
beams [9ll581[59]. In the scattering problem, a probe beam interrogates a (nominally) static liquid, 
with the scattered beam then detected and analyzed. 

Here, we assume that the source of excited hydrodynamic modes within the vortex sheet derives 
from the feature generating sheet-scale vorticity, e.g., a solid or fluid body moving relative to the 
initially static liquid. In contrast to the beam scattering problem [9ll581[59]. we must include the 
space- and time-dependent background velocity field, u = u{y, t)ex, produced by the vortex sheet. 

Model assumptions are as follows: 

a) Due to homogeneity of downstream (x-direction) and cross-stream (z-direction) boundary 
conditions at the moving, vorticity-generating boundary or fluid interface, we assume that 
hydrodynamic fluctuations vary only in the cross-sheet {y-) direction. 

b) For the same reason, we assume that only one hydrodynamic shear momentum current, i.e., 
mass weighted vorticity component, /^^(y,^) = pVxv • e^, appears and that it is directed in 
the cross-stream (2;-) direction. 

c) Consistent with the discussion above, we assume that hydrodynamic fluctuations take place 
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on time and length scales that are short relative to those of the vortex sheet, and long relative 
to do and Td- 

Thus, express the field variables, s, p, w, and ^z, as the superposition of a slowly varying, 
sheet-scale component, f{y,t), and a weak fluctuating sub-sheet component, f'{y,t) : 

f{y,y,t,t) = f{y,i) + f{y,t) (86) 

where s, p, and w = /oV • v = pV • v' are, respectively, the entropy, pressure, and dilatational 
momentum currents. In order to expose sub-sheet-scale processes, sheet-scale time and position 
coordinates, y and i, are magnified using 

y = ey 
i= eit 

where y and t are hydrodynamic-scale coordinates and where 

e = dc/5i 
€i = Tci/{dc/ao) 
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Thus, on the sub-sheet, hydrodynamic scale, time and space derivatives of / are given by: 



dy dy dy 



df df df 



In other words, on the sub-sheet scale, variations in sheet-scale fields are small. 

Extending the approach in [5S] to the problem of hydrodynamic fluctuations within a planar 
vortex sheet, we derive the following system of equations governing s', p', w' , and /x!, : 

[| _ (A; - l)aTV2] p'{y,t) + alw'iy,t) - pp~\k - l)aTV^s'{y,t) = Fi = 0(e,ei) 

[|_(4^ + ^^)V2]^„'(y,t) + VV(y,t)= F2 =0(e,6i) 

(88) 

[d.-aTV^]s'{y,t)-^p-^aTV^p'{y,t)= F3 =0(e,ei) 
[|-z.V2]^;(y,t)= F4 =0(e,6i) 

where k = Cp/Cy is the ratio of specific heats, a-p is the thermal diffusivity, a^ is the adiabatic 
sound speed, /3 is the coefficient of thermal expansion, and v = /i/p and vb = P^bIp-, where p and 
Pb are the shear and bulk viscosity coefficients, respectively. 

Prior to listing the 0(e,ei) terms on the right above, we define 






— f)^ 

-{p'f + T'p)— + A'-7r + A-7v' 



(89) 



where Ajj = (vij + fj,i)/2 is the rate of deformation tensor and ttjj = 2r/iAjj + rj2Aij5ij is the 
stress tensor, and rji and 7^2 are the shear viscosity and dilatational viscosity, respectively (with 
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rji = fj, and 2771+772 = ^^J■ + ^J■B)■ Primes on Aij and vTy denote tensors in the hydrodynamic velocity 
component; here, only derivatives in the cross-sheet (y— ) direction appear. It is seen that products 
of primed and unprimed versions of A and tt represent viscous dissipation due to weak 0(e,ei) 
interaction between hydrodynamic-scale and sheet-scale velocity fields. 

Given Fq, the 0(e,ei) terms above can be expressed as: Fi = {k — l)pf3^^Fo, F2 = — vt • p' — 
pVj^i'Vij, F3 = Fo, and F4 = -p^^Vx (p'vt) - Vx [v' • Vv]. 

Importantly, the 0(1) left hand side of (j88p is, with the exception of the appearance of one 
rather than two shear momentum currents, identical to that describing the hydrodynamic response 
of a nominally static liquid to beam scattering. Thus, focusing on the leading order problem, we 
can immediately apply the results of the scattering problem to sub-sheet organization in unstrained 
vortex sheets. Taking the Laplace-Fourier transform of the above system 



CXD />00 



f'k= I / e''^'^e-^V'(x,t)dxdi 

'0 J~oo 



with x = ye-y, and enforcing the solvability condition leads to 



[s + {k - l)aTk'^] al pp-\k - l)aTk'^ 

-k'^ [s + uik'^] 



^lau2 








[s + axk"^] 











[s + uk'^ 



(90) 



From the original system (or (j90p ). it is clear that the cross-sheet shear momentum current is 
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uncoupled from the other three modes; from ()90p . the associated dispersion relation 

S4 = -uk'^ (91) 

shows that this mode is purely diffusive. The remaining modes follow from (j90p and are given 
by [581E9]: 

si = -ark"^ (92) 

S2,3 = ±iaok - Aok"^ (93) 

where Aq = [vi + {k — l)aT]/2, with vi = ^v + vb- 

Thus, from (j92]) and (f93|) . we surmise that, beyond cross-sheet diffusion of clump-scale vorticity, 
two additional sub-sheet organizing mechanisms act (respectively) : thermal smoothing via a purely 
diffusive cross-sheet entropy mode, and acoustic smoothing via weakly damped, oppositely-directed, 
cross-sheet acoustic modes (with the rate of damping determined by Aq). 

Physically, in the cross-sheet direction, acoustic dilatational waves tend to push out-of-plane 
vorticity fluctuations into a planarized configuration; likewise, within the in-sheet plane, wave re- 
flection between neighboring vorticity fluctuations tend to align the fluctuations in the z — direction. 
Superposed on these acoustic mechanisms, diffusional smoothing of clump-scale vorticity as well 
diffusion of friction-generated thermal energy provide additional, strongly organizing effects. The 
mechanisms are depicted in figure 5; see the caption for further discussion. 
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Figure 5: Organization of disordered clump-scale vorticity is mediated by three hydrodynamic 
modes: a diffusive, cross-sheet entropy (thermal) mode (upper left), a diffusive, cross-sheet shear 
(vorticity) mode (upper right), and a pair of counter-propagating, weakly damped, cross-sheet 
acoustic modes (bottom) . We depict the thermal mode in terms of the thermal degradation (dis- 
persion) of in-sheet clumps. The vorticity mode can likewise be viewed as the hydrodynamic- 
time-scale clump-to-clump transfer of rotational inertia. On these time scales, clumps experience 
significant dispersion so that rotational transport takes place by both the cog-like interaction of 
adjacent 'thermal bumps', and gross intermixing of neighboring clumps. 

6.4.2. Elemental vortex sheets defined 

Having developed a physical picture of sub-BVS vorticity, we can identify a reasonable physical 
embodiment of the sub-BVS scale stochastic processes defined by equation (j79p : elemental vortex 
sheets are defined as planar vorticity layers of thickness 0{6ci), which evolve stochastically within 
any given BVS, on the long BVS time-scale. 

Two features suggest that EVS's can be taken as real, rather than quasi-physical, entities. First, 
since hydrodynamic acoustic, shear, and entropy modes organize, on long BVS length and time 
scales, disordered small- (clump-) scale vorticity, it is clear that these long-acting modes are likewise 
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capable of organizing 0{5c) layers of vorticity. Thus, on BVS time scales, it is physically reasonable 
to view, at any instant, a Burger's vortex sheet, or any vortex sheet, as being composed of 0{6i/6ci) 
elemental vortex sheets. Additionally, it is likewise reasonable to expect that due to long-time-scale 
organization of short-time-scale vorticity, 0{6ci) layers of vorticity evolve stochastically, en mass, 
over long BVS time increments. At. 

6.5. Vortex sheet evolution in non-random strain rate fields 

Prior to discussing vortex sheets in random strain rate fields, we briefly consider physically inter- 
esting aspects of the solution which apply to non-random strain fields. Thus, using (|77|) in (|66p . 
one can write down the vorticity field produced by N^ initially discrete Burger's vortex sheets: 

n(x,t)=y exp — ^ — —- — ^, (94) 

where again, Avi is the incremental change in the streamwise velocity at Xj. 

Several observations follow from this solution. First, from (j80|) and (j8ip . the mean position and 
diffusive spread of the i sheet, viewed as the stochastic evolution of a cloud of elemental vortex 
sheets, are given respectively by: 

< Xi{t) >x= Xi exp {-h{t)) (95) 

and 

< {x^{t)- < Xi{t) >x? >x= 2up{t) (96) 
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where h{t) and p{t) are given by (j72p and (j76p . Thus, in non-random strain rate fields (and due to 
the hnearity of the problem), each sheet (in the collection of A'^j,' sheets) migrates toward x = 0, each 
diffusively spreading at a spatially uniform rate; at large times, all N^ sheets coalesce at x = 0, with 



the thickness of the collection, y^2h'p(t), identical to that of a single sheet. In the case where k{t) is 
constant, the above mean and variance simplify to those for the equivalent OU process considered 



in [8]. Likewise, identifying a/< (Xii't)— < X-ii't) >a;)^ >x as the thickness of sheet i, we obtain the 
constant strain rate expressions for sheet position and thickness given, e.g., in [52y53j. 

Second, and considering the long-time strength of the composite sheet in the case where k is 
either (a positive) constant or approaches a long-time limit, k — )• koo, h{t) — )• oo and p{t) — )• l/{2koo), 
and the strength of the coalesced collection of Nx sheets, from ()94p . is the sum of the initial 
strengths: 

Thus, for example, when the strength of one of the initial sheets dominates all others, corresponding 
to a locally large change in the initial stream-wise velocity profile, (|97p shows that the final sheet 
strength is largely determined by the strength of this sheet. In contrast, in situations where the 
initial streamwise velocity profile takes on positive and negative magnitudes, and more specifically, 
is characterized by a limited number of Fourier modes, (197p shows that due to vorticity cancellation 
the final strength is small. 
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6.6. Single and continuous sets of Burger's vortex sheets in random strain rate 
fields 

Focusing first on the dynamics of a single vortex slieet (N^ = 1 in (j94p ). we note tliat since 
equations ([M|) . (j95|) . and ([96j) apply to individual realizations of k{t), then the corresponding 
ensemble averaged instantaneous vorticity, sheet position, and sheet thickness, taken over k, can 
be expressed respectively as: 



(Ox, t)) = ( , exp \^—, — V- — -] ) (98) 

^ ^ ^ ^^^up{t) ^ 47ri/p(t) J/ ^ ^ 



(<x(i)>x> = a;o<exp(-/i(i))> (99) 

and 

( < {X{t)- < Xit) >.f >x > = 2i^{p{t)) (100) 

where h{t) and p{t) are again given by ()72p and ()76p . and Af is the incremental (step-function) 
change in the streamwise velocity at Xq- 

Explicit formulas for these quantities can be obtained when the strain rate field is statistically 
stationary, has positive mean, and a gaussian random component: 

k{t) = ko + k'{t) (101) 

where, for stability, we assume a condition like k^ > a({k')'^\ where a is an undetermined positive 
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number. Under these conditions [60] 



exp[- [ k'{s)ds])=exp[l- [ [ {k'{s)k'{s'))dsds'] (102) 

Jo ' 2 Jo Jo 



Inserting the correlation function, i?, and using the property 



i?(|r|) = (A;'(i' + r)A:'(t')> (103) 



appropriate to stationary processes, allows restatement of (|102p in the form: 



exp[- [ k'{s)ds]) =exp[l [ [ R{\s - s'\)dsds'] (104) 

^0 ' 2 Jo Jo 



In turn, the identity [53 



Iff R{\s - s'\)dsds' = f {t- s)R{\s\)ds (105) 

2 Jo Jo ^0 

is useful when evaluating the right side of ()104p . 

Thus, inserting (|104p and (|105p in (|99p yields the following formula for the (single) sheet's 
ensemble average position, applicable to random strain fields of the form in (jlOip : 



< Xit) >x ) = Xoe-^°'exp [ [ {t- s)R{\s\)ds] 

Jo 



(106) 
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In order to obtain the corresponding single sheet spread, we use ([76]) to write 



(p(t)) = e-2'=°*(exp [ - 2 /" k'{t')dt'] ■ I (e^^^"*' exp [2 / k'{s')ds']^dt') (107) 



and assume that the two terms separated by the center dot represent independent random variables: 



/■t /"t rt 

{p{t)) = e-^'''''{exp[-2 k\t')dt'])-{ (e^'^"*' exp [2 / k'{s')ds']yt') (108) 



or 

ft ft , i-t' 



(p(t)) = e-^^oti^ exp [ - 2 / fc'(t')dt'] > • / (e''^"*' ( exp [2 / }^{s')ds'\ >) dt' (109) 



JO ^ JO 



Using p04|) and (|105|) to evaluate the expectations in (|109p and inserting the results in (|100|) then 
yields: 

( < (X(i)- < X{t) >.f >. > = 21.6-2'=°* exp [4y" {t-s)R{\s\)ds]- j' (e^^"*' exp [4^* (t'-s')^(l^'l)ds'])rft' 

(110) 
Considering, for illustration, the case where the random strain field is delta correlated in time 

{k'{t' + T)k'{t')) = k6{T) (111) 

where /c is a positive constant (determined, e.g., by experiment), one finds that sheet spread no 
longer approaches a fixed long-time limit (as in the case of deterministic strain rate), but rather 
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exhibits exponential, i.e., superdiffusive [6T] growth: 

( < {X{t)- < X{t) >.? >x ) = TT^exp (m[l - exp ( - {2K + 4fe)i)]) (112) 

In the limit of weak random strain, specifically kt « 1 and ko » k, (]112p simplifies to the 
constant strain (non-random) solution [8]: 

< < ixit)- < x{t) >xf >0 = f [1 - exp ( - 2kot)] (113) 

Kg 

[See section 16.7.11 below for a brief description of the conditions necessary for introducing an as- 
sumption of delta correlated strain rates.] 

6.6.1. Continuous collections of Burger's vortex sheets: Feynman-Kac solution 

Turning to the case of multiple Burger's vortex sheets, we consider the general case where the initial 
streamwise velocity component exhibits a continuous variation in the x— direction, 

v{x,0) = Vo{x) (114) 

Here, a representative solution for the mean continuous vorticity, (a;(3;,i)), can be obtained using 
the Feynman-Kac solution [20] of the backward-time version of the vorticity equation (j62p : 

{n{x,t)) = (i^^(i)=,,,[f](x(0),0)]> • (e'^W) (115) 
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where -E'^(t),t denotes the expectation associated with the stochastic process x samphng the initial 
condition, O(x,0). 

As a quick check of this formula, we note that when k{t) = ko is non-random and constant and 
Vo = CqX, i.e., the initial streamwise velocity is of Couette form, the solution in pi5|) yields 

17(x, t) = Coe^°* (116) 



This solution in turn satisfies the vorticity transport equation, (j62p . along with the initial condition, 
i7(x, 0) = Co- Equation (J116p holds under conditions where the time scale for setting up the Couette 
profile, say r, is short relative to the strain field time scale, r^ = k^^. On the longer scale, diffusion of 
vorticity saturates out quickly and the appropriate initial condition, to order k^^ /t, is ri(x, 0) = Cq. 
Under these conditions, the initial vorticity, Cq, is exponentially amplified by stretching. 

As shown immediately below, in the case where k{t) is random and Vo again equals CqX, the 
Feynman-Kac solution in (|115|) coincides with the leading order low- viscosity (high Reynolds num- 
ber) forward-time solution of (j62p . The paradox of the viscous solution coinciding with the inviscid 
solution is resolved by requiring (in analogy with the deterministic strain case immediately above), 
that for every realization of k(t), the initial (linear) streamwise velocity profile forms on time scales 
much shorter than say, \k~^ {t)\max] for each realization, diffusive effects again saturate out over the 
short time scale. We now briefly consider the latter, low viscosity problem. 
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6.7. Single-sheet and collective behavior under inviscid conditions 

In the case of single sheets, significant questions include the effects of random strain on time- 
varying mean (ensemble averaged) sheet strength, position, and spread, while in the case of sheet 
collections, the key question centers on evolution of the ensemble average (position and time- 
dependent) vorticity. Attention is again focused on the case where k{t) is given by (jlOip with k'{t) 
again a stationary zero-mean gaussian process. 

Considering first single sheets, the position of the sheet under the action of an individual real- 
ization of k{t), 

Xit) = e-^Wx(O) (117) 

follows either from the SDE governing the motion of elemental vortex sheets, (j79p . with the stochas- 
tic term suppressed, or by noting from the MOC solution for Q{x,t) (using the inviscid form of 
(j62p ) that any given sheet evolves along characteristics given by (|70p . In the case of elemental 
sheets, since no diffusion takes place, all EVS's within the Burger's sheet evolve identically under 
the action of any given realization of k(t). In the following, we express the sheet's initial position 
as x(0) = Xo- 

Given x(0) the ensemble average sheet position over the random strain rate field can be calcu- 
lated using (fTOil) and (fT05]l : 



(x(t))=Xoe-'^°*<exp[- / k'{s)ds]) 

Jo 



XoC ''°*exp 



{t - s)R{\s\)ds 





:ii8) 



Since (i — s)i?(|s|) is positive for s < t, then it is seen that any stationary random strain rate field 



56 



inhibits sheet migration toward the origin, x = 0. 

Given (^xi't)), the sheet spread (or squared thickness), o"^, can again be calculated: 

4 = (Mt) - {x)f) = xle-''^' [{ exp -2 j^ k\t')dt') - ( exp - f^ k'{t')dt'f] 



(119) 



or 



al = xle-^^'''\ 



X o^ 



[exp[4 / R{\s\){t - s)ds] - ex.p[2 R{\s\){t - s)ds]\ (120) 



6.7.1. Inviscid vortex evolution in a rapidly decorrellating random strain rate field 

This example briefly considers the spread of a single Burger's vortex sheet under conditions where 
the correlation time, Tc = L {k' {t)k' (0)) dt/var{k'), of the random strain field, k'{t), is short 
relative to the time scale, Tb, associated with the vorticity-inducing flow, v{x,t). Thus, deflne a 
normalized correlation function ri(T) = {k'{T)k'{0)) /cf, where cf is a normalization factor obtained 
from J_ rj{T)dT = 1. In addition, let e = Tc/ti, << 1, and renormalize the correlation function to 
the longer BVS time scale: ??e(f) = e^'^rj{e^^T), where f = r/e. Since r/e(f) — )• (5(f) as e — )• (and 
f = 0{Tb)), then it is clear that the assumption of delta correlated statistics on the BVS time scale 
becomes increasingly valid as the ratio of the Tc/rb becomes increasingly small. 
Thus, let k'{t) be delta correlated. In this case, (jllSp yields: 

{xit)) = Xoe(-'^°+^)* (121) 

Although this suggests that for large enough random strain, k > ko, the mean position of the BVS 
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can track away from the origin, x = 0, this cannot be confirmed without detailed consideration of 
sheet stabihty under random strain rates. 

Under the same assumption, a variety of interesting long-time behaviors emerge, depending on 
the size of k relative to ko '■ 

al = xle-^^'''[e^~''-e^'^'] (122) 

In particular: 

i) when k < ko/2, corresponding to mean strain dominating the random component, sheet 
thickness goes to (i.e., o"^ — )■ 0) as t — ;• oo; 

ii) when k > ko/2, sheet thickness grows without bound (cj^ — )■ oo) as t — )■ oo; 

iii) in the special case where k = ko/2, a fixed, non-zero thickness can be achieved, a^. — t- Xo as 
t — 7- oo, indicating, both here and in the other cases, that the random strain functions as an 
effective agent for diffusion. 

In each case, the increase in sheet spread amplitude with initial position, x^, refiects the increase 
in random velocity amplitude, v'{x,t) = k'{t)x', with x. Finally, and once again, cases ii) and iii) 
are subject to the proviso concerning sheet stability in increasingly random strain fields. 
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6.7.2. Single sheet and continuous vorticity in the inviscid limit 

Considering the averaged evolution of the vorticity, (^Q{x,t)'j, we first introduce the limit z^ — t- 
into the single sheet solution (j77p (which is applicable to individual realizations of k{t)), to obtain 

n{x,t) = Av6{x - Xoe-^(*)) (123) 

where again h{t) is in defined in (|72p . Physically, a vortex sheet initiated at Xq simply advects with 
the random strain field. Expressing the delta function in its Fourier representation form, followed 
by evaluation of the ensemble average then yields: 



{n{x,t)) = — e'"'''(exp[-imxoe-''°^e-^o'''(^')'i^')dm (124) 

27r J_^ 



We refrain from analyzing this expression in depth, but note in the case where mean strain 
dominates random strain, kg >> k, the long-time limit of (|124p yields 



[nix,t)) = / e^^^dm = Av6ix) (125) 

27r ./_^ 



Thus, when strain rate is strongly deterministic, and under inviscid conditions, the initial, infinitely 
thin vortex sheet remains infinitely thin, migrates to x = 0, and becomes infinitely strong, all 
physically reasonable results. 

Turning next to the inviscid evolution of a continuous set of vortex sheets, and focusing initially 
on a single realization of random k{t), we first write the discrete multi-sheet solution corresponding 
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to dMl) as 

n{x, t) = Y^ Avi6{x - Xie~^(*)) (126) 



i=l 



and then re-express this as an integral: 



n{x,t) = r ^}^^li^6ix - x'e'^^'^)dx' (127) 



oo '^•^ 



Integrating then gives 



O(x,i) = e^W^t;(x,™t(t),0) (128) 

ox 



where 

Xinit{t) = xe^(*) (129) 

is the initial x position (at t = 0) of the characteristic passing through x at time f; refer to equation 

()117|) . Thus, for any given realization of k, the initial vorticity, generated at Xinuit) and having 

magnitude 

d 

is amphfied (via stretching) by a factor of e^*'' . 
The ensemble average vorticity assumes the form 



d_ 
dx' 



{n{x,t)) = (e'^W— t;(x.™t(t),0)> (130) 



The consistency of this inviscid solution is shown by noting that the representative general (viscous) 
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solution obtained via the Feynman-Kac approach, equation (jllSp . simplifies to (jl30p when the 
stochastic (i.e., viscous) term in the SDE governing x(i) is set to zero. In other words, for any 
realization of k{t), the stochastic processes used to construct the FK solution track along the inviscid 
characteristics, dx = —k{t)xdt, so that 

d 

yielding (fTSOll . 

As a simple example, in the case where the initial streamwise velocity is of Couette form, 
v{x',t = 0) = i^ox' , where Oo(= Cq above) is constant: 

(n(x,t)) =J7o(e'^W> (131) 

When k{t) is given by (jlOip with, e.g., a delta correlated random component k'{t), (jlSip assumes 
the form: 

(!^(x,t)> =J7oe(^°+^)* (132) 

Hence, in this simple example, under inviscid (i.e., high Reynolds number) conditions, continuous, 
initially constant vorticity is amplified not only by the mean strain field, ko, but also by the random 
component. 
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7. Concluding remarks 

Three themes are pursued. First, by comparing the representative stochastic solution of a hnear, 
nonhomogeneous, drift-diffusion problem against the corresponding Green's function solution, we 
obtain a set of equalities relating stochastic expectations in the former to Green's function convo- 
lutions in the latter. Importantly, these equalities expose a framework within which stochastic and 
Green's function methods can be applied in concert to a range of problems. In broad terms, and 
as illustrated above, the framework allows exploitation of two generic constructs - delta functions 
and Wiener processes - for mathematically and physically modeling and probing an array of linear 
and nonlinear, deterministic and random problems. 

The second theme centers on testing application of time- incremental GF's to solution of nonlinear 
evolution problems. Two canonical problems, Burger's equation and the nonlinear Schrodinger 
equation, provide test beds. The first Test Case, focused on the simplest embodiment of nonlinear 
drift-diffusion problems, leads to the following observations: 

i) Transforming from an incremental to nonincremental solution rests on elimination of nonlin- 
earity from the incremental GF. Here, the transformation is accomplished via the following 
procedure: 

(a) guess a transform, 0(m), parameterized in the unknown u, 

(b) force the evolution of (/>(u) to be purely diffusive, and 

(c) attempt to match the evolution of u that emerges from the diffusion of 4>{u), equation 
(fSHI) . to the original nonlinear evolution equation governing u, equation (fT9]) . 
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ii) Proper choice of the time step, embodied, for example, in equations (j27p and ()28p . allows 
approximate solution of the otherwise difficult-to-solve backward adjoint problem. This step 
is essential to obtaining an analytical incremental GF. 

iii) The procedure outlined in i) and ii) allows derivation of the Cole-Hopf ansatz. Moreover, we 
anticipate that a similar approach can applied to other nonlinear evolution problems. 

The second Test Case, an example of nonlinear parabolic wave propagation, illustrates application 
of asymptotic approaches for transforming from incremental to non-incremental solutions; here, a 
well-known soliton solution is recovered. 

The last theme revolves around use of the GFSM as a tool for probing physical features in 
problems characterized by some element of randomness. Here, a thorough investigation of single, 
multiple, and continuous sets of Burger's vortex sheets evolving in deterministic and random strain 
rate fields, under viscous and inviscid flow conditions, is presented. The main results are as follows: 

i) For delta function initial conditions, the vorticity transport equation assumes the form of a 
Fokker-Planck equation of an Ornstein-Uhlenbeck process. 

ii) Correspondingly, the evolution of a Burger's vortex sheet, i.e., the movement of its mean 
position and its time- varying spread, can be viewed as the evolution of a constituent cloud of 
elemental vortex sheets, governed by an OU stochastic differential equation. 

iii) A physical picture of EVS's is developed by first focusing on the short time-scale vorticity 
evolution of particle clumps. On short acoustic time scales, sub-sheet vorticity is three- 
dimensional, disordered, and strongly diffusive. 
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iv) Considering the fundamental question of how disordered, clump-scale vorticity becomes or- 
ganized on long BVS time scales, we study an analog problem of organization within non- 
strained planar vortex sheets. In this case, an analysis of sub-sheet hydrodynamic modes 
suggests that organization is mediated by a combination of damped, cross-sheet acoustic 
modes, a diffusional cross-sheet shear mode, and a diffusional cross-sheet entropy mode. 

v) A number of analytical results describing the motion and spread of individual, discrete col- 
lections, and continuous sets of Burger's vortex sheets, evolving within deterministic and 
random strain rate fields, under both viscous and inviscid conditions, are also obtained. 

Appendix: Definitions; notes regarding incremental GF's 

Definitions of forward and backward time coordinates and forward and backward evolution and ad- 
joint equations are first given. Four important points regarding determination of time-incremental 
Green's functions are then highlighted. 



Forward and back^vard time coordinates and equations 

Coordinates in the forward, i.e., physical time direction are denoted as t or t'; the backward time 
coordinate is denoted as s or s', A drift-diffusion (or drift-diffusion-like) evolution problem is of 
forward time form if: i) the signs on the time derivative and diffusion terms in equation ([TJ differ, 
and ii) the evolution is initiated from some known initial state. Likewise, an evolution problem is 
of backward time form if: i) signs on the time derivative and diffusion terms in ([T|) are the same, 
and ii) the evolution ends at some known final state. 
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Notes on incremental GF's 

We highlight four points. 

First, the incremental Green's function, G, must meet the following requirements: 
i) It must satisfy the adjoint problem associated with the evolution problem of interest, 
ii) For forward form adjoint problems, again associated with backward time evolution problems 
(both stated in terms of the backward time coordinate, s'), G must behave as 6{x.' — x) as s' G 
[sj,Sj+i) approaches Sj; here, x' and x are variable and fixed points. Likewise, for backward form 
adjoint problems, associated with forward time evolution problems, both stated in terms of the 
forward time coordinate, t' , G — )• J(x' — x) as t' — )• tj (where tj < t' < ij+i). Note that, e.g.. 
Barton |29] provides a useful compendium of delta function representations. 

iii) In problems where boundary effects are important, G must, over each time interval, satisfy 
appropriate homogeneous Dirichlet or Neumann (or possibly) mixed boundary conditions. 

Second, from ([7]), it is apparent that the adjoint equation associated with the forward time 
evolution problem in (l2])-(j3]) is of backward time form in the forward time coordinate, t' (i.e., 
the signs on the diffusion and time derivative terms are of the same sign). Likewise, adjoint 
equations associated with backward evolution equations are of forward form in the backward time 
coordinate s'. Importantly, while it is often difficult to determine analytical solutions to backward 
time problems |361I37] . we can sometimes exploit the smallness of At' or As' to adapt relatively 
simple forward time solutions to the approximate solution of associated backward problems. This 
point is illustrated in our solution of Burger's equation in Test Case 1. 

Third, a variety of methods are available for computing non-incremental Green's functions 
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|28p29j. all of which can be applied to computing incremental solutions. Barton [29] provides 
an accessible description of a number of approaches, including, e.g., eigenfunction expansions for 
problems subject to boundary effects. We note in passing that numerical application of incremental 
Green's function solutions has antecedents in the literature on the boundary element method [M] 
and the so-called Green element method |18y35j. The present framework (GFSM), however, differs 
from these in an essential way due to its combined reliance on Green's function and stochastic 
process methods. 

Fourth, when interpreting incremental Green's functions as incremental transition densities, 
incremental G should, strictly speaking, satisfy two consistency conditions: it must allow recovery 
of the local mean drift, b(x , s ), and the local diffusion matrix, B(x , s ) 

b(x , s ) = lim -— - / (x — X )G(x , s x , s + As )ax (A-1) 

As'-^O As' Jo 



and 



Bij{-K ,s] 



lim -— - / (Xj —xA[x, — x)G(x,s|x ,s + As )ax , (A-2) 

As'^o As' Jd ■> ■> 

where in this paper, B^j = 2i'5ij, and where |x — x | < e, with e arbitrarily small [8]. 
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